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Abstract 

Grazing-incidence X-ray diffraction data are combined with energy-minimization calcula- 
tions to analyse the atomic structure of the Al-rich (3\/3 x 3\/3)R30° reconstructed surface 
of sapphire a— Al2O3(0001). The experiments on the BM32 beamline of the ESRF provide 
the non-integer-order diffraction intensities and, after Fourier transform, an incomplete Pat- 
terson map. The computer simulations are implemented to obtain structural information 
from this map. In the simulations, the interactions between the Al overlayer atoms were 
described with the Sutton-Chen potential and the interactions between the overlayer and 
the sapphire substrate with a laterally modulated Lennard-Jones potential. We have shown 
that the hexagonal reconstructed unit cell is composed of triangles where the two layers of 
Al adatoms are FCC(lll) ordered whereas between the triangles the stacking is FCC(OOl). 
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1 Introduction 

The (0001) surface of sapphire a— AI2O3 is, because of its insulating character, an 
important substrate for very high frequency microelectronic devices. Therefore it 
is of crucial importance to understand the stability and structure of this surface at 
high temperatures. In equilibrium at room temperature, a clean (0001) surface is 
unreconstructed and is terminated with a single Al layer on top of the last oxygen 
layer [1-3]. Upon heating up to ~1500 C in ultra-high vacuum (UHV), oxygen 
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evaporates from the crystal, leaving an Al rich and reconstructed surface behind it 
[4]. Al atoms that are left on the sapphire surface are responsible for a sequence of 
surface reconstructions: (\/3 x v^)R30° when sapphire is heated to about 1 100 C, 
(2^3 X 2V3)R30° at 1150 C, {3^3 x 3y3)R30° at 1250 C and finally (\/3T x 
a/3T)R±9° structure at 1350 C [5]. So far only the structure of the most stable, i.e., 
of the {\/31 X \/3T)R±9° reconstruction has been satisfactorily determined [6,7]. 
Very recently, this structure has also been investigated by Barth and Reichling with 
the atomic-scale-resolution scanning force microscope [8]. No information on the 
real-space structures of other reconstructions is available so far. 

In this letter we concentrate on the {sVSx 3\/3)R30° reconstruction of the sapphire 
(0001) surface. Although this is a transient reconstruction, it was prepared in a very 
well defined state with large domain sizes [5,9]. 

Our analysis is based on grazing-incidence X-ray diffraction (GIXD) data where 
some surface diffraction intensities are hidden in the bulk Bragg peaks and the 
corresponding Patterson map is incomplete [10]. To obtain the real-space structure, 
therefore, we combined the Patterson map with computer simulations using semi- 
empirical potentials. 



2 Experiment 

The structure was investigated experimentally by grazing-incidence X-ray diffrac- 
tion using the SUV set-up [1 1] of the BM32 beamline at the ESRF. This is equipped 
with a large UHV chamber (3 x 10^^^ mbar base pressure) with RHEED, AES, 
mass analyser, a high temperature (up to 1500 C) furnace and several deposition 
sources. The sample was prepared by heating at 1250 C for 20 min, which yielded 
a (3v^ x 3v^) reconstruction of very high quality, as checked by in situ RHEED, 
and negligible segregation of contaminants, as checked by AES. A 18 keV doubly 
focused monochromatic X-ray beam was used, with a 0.2(H) x 0.2(V) mm^ size 
and 0.5(H) x O.l(V) mrad^ convergence at the vertical sample location. The an- 
gular acceptance of the Nal detector was set to 20(H) x 4(V) mrad^. The incident 
angle was fixed at the critical angle for total external reflection during the whole 
measurement. A total of 1790 in-plane reflections were measured, showing p6mm 
symmetry. These were integrated over azimuthal rotation and corrected for back- 
ground, effective area, polarisation and Lorentz factors [12] to extract the in-plane 
structure-factor amplitudes Fhko- This reduced to 240 inequivalent peaks. The ex- 
perimental diffraction pattern in 1/6 of the reciprocal space is shown as right-hand 
semi-circles in Fig. 1, with their radii proportional to Fhko- The h and k indices are 
in reciprocal lattice units of the (3a/3 x 3 v/3)R30° surface reconstruction. Since the 
phases of F^ko are unknown, we cannot directly determine the real-space density 
distribution. In addition, the amplitudes at the bulk Bragg reflections are unknown 
since the contribution of the surface reconstruction is hidden in the bulk contribu- 
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Fig. 1. Diffraction pattern of the (3-\/3 x 3-\/3)R30° reconstructed surface of sapphire 
(0001). The areas of semi-circles are proportional to the intensity of diffraction spots. The 
right-hand semi-circles are experimental results and the left-hand ones correspond to the 
structure proposed in this paper. The in-plane bulk reflections are omitted from the graph. 
The dashed lines show the unreconstructed reciprocal cell. 

tion. Therefore the Fourier transform of the experimental F^^^q gives an incomplete 
Patterson map, i.e., an incomplete electron (atom) density-density correlation func- 
tion. The incomplete Patterson map of the (3v^ x 3 v^)R30° reconstructed sapphire 
(0001) is shown in Fig. 2(a). 



3 Simulations 



Although the Patterson map does not give the real-space structure, it is the starting 
point to build the real-space models. It tells us that the surface atoms that contribute 
to the more compact regions of the Patterson map are ordered on a hexagonal lat- 
tice, very similar to the "ordered domains" seen also in the {^/Sl x \/3T)R±9° 
reconstruction of sapphire [6,7]. Therefore we constructed the overlay er assuming 
that sapphire is covered by two FCC(lll) layers of Al with the nearest-neighbour 
separation ~ 4 % larger than the corresponding sapphire spacing. The lattice misfit 
inevitably leads to perturbations in the close-packed FCC(1 11) structure, creating a 
commensurate structure composed of "domains", where the overlayer atoms have 
perfect FCC(1 1 1) or hexagonal close-packed ordering, separated by a kind of "do- 
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Fig. 2. The Patterson map of one reconstructed unit cell on sapphire, (a) Incomplete Pat- 
terson map based on the experimental diffraction data only, (b) Patterson map when the 
intensities calculated from the model at the in-plane bulk reflections are added to the ex- 
perimental diffraction data. 

main walls". Whereas the atomic structure in the domain walls is disordered in the 
{Vsi X Vsi) reconstruction [6-8], the Patterson map of the (3^3 X 3^3) recon- 
struction in Fig. 2 suggests a much more ordered structure of domain walls. We 
generated several initial structures with hexagonal close-packed atom ordering in 
the domains with domain walls between them. 

Such reconstructions are too large to be investigated with ab initio methods. There- 
fore we described the interactions between the Al overlayer atoms with the semi- 
empirical Sutton-Chen potential [13,14] and the effect of the substrate on the over- 
layer with a phenomenological potential field with hexagonal symmetry. The en- 
ergy minimization was described in detail in [7]. 

The semi-empirical many-body Sutton-Chen potential 



includes the core repulsion potential (first term) and the bonding energy mediated 
by the electrons (second term). Vij is the separation between the atoms i and j, pi 
is an effective local electron density at the site i: 



and a is the lattice constant of an Al FCC crystal, e and C are parameters of the 
model which, together with the exponents n and m, determine the repulsive and 
cohesive energies, respectively. We used the following values for the potential pa- 
rameters of Al [13]: 

m = 6, n = 7, e = 33.147 meV, C = 16.399. 
The potential was truncated continuously (with a fifth order polynomial) between 
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r/ro = 3.17 and 3.32 (ro is the nearest-neighbour distance). In this way the inter- 
action with 68 neighbours would be included if Al were perfectly ordered in two 
FCC( 111) planes. 

In general, the overlayers or adatoms are weakly bound to the substrate and well 

separated from it [15]. In addition, the substrate is much stiffer than the Al over- 
layer, therefore the relaxation of the substrate caused by the overlayer was neglected 
in the simulations. The substrate potential was expanded in a power series and only 
the six lowest-order terms were retained [7], 



Us = Ulj{z) 



Ul — cos{ki ■ f) cos(/c2 • r) cos(A;3 
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are the unit vectors in the plane of the surface and ag the substrate lattice constant. 
The parameter U l controls the strength of the lateral modulation of the potential, 

and Ulj, 



Ulj{z) = Uo 



- -2 - 



zo 



(5) 



its dependence, in the direction perpendicular to the surface. Expression (5) is 
the Lennard- Jones potential integrated over the semi-infinite substrate where Uq is 
the depth of the substrate potential and zq determines the position of the minimum 
and the width of the potential in the vertical direction. The potential has six equally 
deep minima in the surface unit cell, to accommodate three overlayer atoms in each 
monolayer. Thus, the overlayer-substrate potential is described by three variational 
parameters: Ul, Uq and Zq. 

The parameters Ul,Uq and zq were varied to minimize the difference between the 
experimental and calculated structure factors [5]: 

2 _ 1 (l-^fefeol ~ l-^i^fcO l) 

hk K^hko) 



where the summation is over N measured diffraction peaks, ctJ^^q are the experi- 
mental errors in Ff^n whereas &n are the structure factors of the trial structures. 
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The initial structures with the number of atoms in reconstructed unit cell varying 
between 134 and 140 were first tested for their stability in a reasonable range of 
substrate potentials and then the parameters Ul, Uq and zq were varied until the 
minimal was reached. In this way we narrowed the possible initial structures 
to one, having 136 surface atoms in the reconstructed unit cell, shown in Fig. 3. 
This structure is very stable, it does not vary after random atomic displacements 
up to 0.286 A and subsequent energy minimization. The variational parameters in 
the minimum of (x^ = 2.54) are: Uo = 0.28 eV, Ul = 23 and Zo = 5.10 A. 
The average Sutton-Chen potential energy per adatom is Esc = —3.045 eV and 
the average adatom-substrate energy Eg — —0.230 eV, giving E — —3.275 eV for 
the average total potential energy per adatom. varies by less than 10% for Uo 
varying between 0.16 and 0.70 eV. 




Fig. 3. The (3\/3 x 3\/3)R30° reconstructed unit cell with 136 atoms (73 in the lower and 
63 in the upper plane) in the calculated model. The darker spheres represent the atoms in 
the lower plane. The underlying triangular net represents the periodicity of the simulated 
substrate potential which has minima in the centres of the triangles. 



4 Results and Discussion 

We would like to emphasize the strength of the present combined experimental- 
theoretical approach compared to the standard surface X-ray diffraction analysis 
based only on the Patterson-map analysis and minimization [16]. 

The first advantage is related to the weakness of the standard method where the 
surface structure factors corresponding to the locations of the bulk contributions 
are not included in the Patterson analysis, thus leading to an incomplete Patterson 
map [10]. Omission of some peaks from the diffraction data can introduce errors 
in the final atomic positions. In Fig. 2 it is demonstrated that the peaks correspond- 
ing to (/i, k) = (0, 0) (mod 9) have a dramatic effect: when they are omitted, 
the incomplete Patterson function shows 8 atoms along the reconstruction primi- 
tive vectors [Fig. 2(a)] instead of 9 that are retrieved if we complete the diffraction 
data with the model values [Fig. 2(b)]. In case of free FCC(lll) stacked Al the 
lattice parameter is ~ 4 % larger than the corresponding sapphire distance. In the 
present case, the misfit is ~ 2 — 3 %, showing that the Al layer is strained by 
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Fig. 4. The relaxed atomic configurations of the (3\/3 x 3\/3)R30° reconstruction. Clearly 
seen are the triangular domains with hexagonal [FCC(lll)] structure separated by "walls" 
where the structure is FCC(OOl). At the comers of the triangles the density is lower. 

the underlying sapphire lattice. Such overlayer yields significant intensities at the 
bulk peaks locations (corresponding to no expansion) as well as in the neighbour- 
ing peaks (12.5 % expansion). A standard analysis would thus favour the second 
components in the Al atom positions, hence probably over-estimating the expan- 
sion of the Al overlayer. The present approach avoids this effect because the bulk 
component is re-introduced via the interaction with the substrate, and thus the final 
structure would be closer to the real one than with the standard analysis. 

The second advantage of the present method is that it yields a high-symmetry struc- 
ture with negligible disorder which, we believe, is close to the "ideal" structure of 
the (3\/3 X 3\/3) reconstruction while the experimental data correspond to a struc- 
ture which may be substantially disordered. 

The final, relaxed structure is shown in Fig. 4 and in Fig. 1 we compare the cal- 
culated (left-hand semi-circles) diffraction pattern with the experiment (right-hand 
semi-circles). Specific for this structure is a triangular pattern with the period SVSag. 
Inside these large triangles, the two Al planes are ordered hexagonally in FCC(1 1 1) 
planes. In the "down" triangle (upper left-hand side of Fig. 3) the lower-layer Al 
atoms are located in the centres of the underlying small up-triangles, represent- 
ing the substrate potential, and the upper-layer atoms in the centres of the small 
down-triangles. In the "up" triangle (lower right-hand side of Fig. 3) the stack- 
ing is reversed. The two large triangles together, thus, form the rhombic recon- 
structed unit cell. Between the triangles, the lower-layer Al atoms are ordered in 
squares; together with the atoms in the upper layer they form semi-octahedra with 
an FCC(OOl) orientation. The semi-octahedral ordering is needed to compensate 
for the ~ 4% lattice misfit between the sapphire substrate and the Al overlayer. In 
the comers of the large triangles we find small hexagonal vacancy islands where 
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several Al atoms are missing in the upper layer. The average density of the over- 
layer atoms is 136/27 ^ 5 Al atoms per unreconstructed surface unit cell and is 
surprisingly close to the density proposed for the (\/3T x \/3T) reconstruction [7]. 
In the large triangles, the nearest-neighbour Al-Al spacing of the lower layer is be- 
tween 2.82 and 2.84 A, giving a contraction of up to 1.5 % with respect to bulk 
Al. The spacing in the FCC(OOl) octahedra is similar, with a contraction of up to 
2 %. In principle, our ground- state simulations predict a chiral structure with p6 
symmetry. The p6mm symmetry is restored either by superimposing both chiral 
orientations or by room-temperature thermal fluctuations since the deviations from 
the mirror symmetry are small, of the order 0.03 A. It is interesting to note that the 
structure of the (3\/3 x 3\/3) reconstruction, as proposed in this paper, is much less 
disordered than the final, (V^ x a/sT) reconstruction [7,8]. 

The difference between the experimental and simulated diffraction patterns in Fig. 
1 might have two origins. First, the adatom-adatom and the adatom-substrate poten- 
tials were treated in an approximate way, and second, there can be some disorder in 
the surfaces. Indeed, experimentally, the reconstruction with 3\/3 symmetry can be 
prepared in very different states of order (as checked by electron or X-ray diffrac- 
tion), depending on the exact annealing temperature and duration. A slight variation 
from the best conditions leads to a high level of background scattering as well as to 
different relative ratios of the structure factors. A minimization of x"^ with respect 
to the Debye-Waller factor using the ROD program [16], for instance, and keeping 
all the atomic coordinates fixed improves the value to = 2.05. The corre- 
sponding Debye-Waller factor increases to cr = 0.25 A, indicating that some kind 
of disorder (i.e., variation of the atomic structure between different reconstruction 
cells) is present on the prepared surface. In contrast, a standard minimization of x'^ 
with the ROD program varying the atomic positions yields a much better agreement 
(X^ as low as 0.7!) but for very disordered and unphysical structures. 

To conclude, on the basis of our GIXD experiments combined with the energy 
minimization using effective potentials, we propose a real-space structure of the 
(3\/3 X 3\/3) reconstruction of Al2O3(0001). We have shown that the combina- 
tion of GIXD experiments with the simulations based on energy-minimization is a 
powerful tool to predict and understand the large unit cell reconstructions, in partic- 
ular on oxides or insulators, where atomic-scale information from scanning force 
microscopy is still extremely scarce. 
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